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We consider the phase coherent transport of a quasi one-dimensional beam of Bose-Einstein 
condensed particles through a disordered potential of length L. Among the possible different types 
of flow we identified [T. Paul, P. Schlagheck, P. Leboeuf and N. Pavloff, Phys. Rev. Lett. 98, 210602 
(2007)], we focus here on the supersonic stationary regime where Anderson localization exists. We 
generalize the diffusion formalism of Dorokhov-Mello-Pereyra-Kumar to include interaction effects, 
ft is shown that interactions modify the localization length and also introduce a length scale L* for 
the disordered region, above which most of the realizations of the random potential lead to time 
dependent flows. A Fokker-Planck equation for the probability density of the transmission coefficient 
that takes this new effect into account is introduced and solved. The theoretical predictions are 
verified numerically for different types of disordered potentials. Experimental scenarios for observing 
our predictions are discussed. 
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I. INTRODUCTION 

The absence of diffusion of waves in disordered media 
was predicted by Anderson 50 years ago Originally 
proposed in the context of electronic transport in disor- 
dered crystals, it has since been observed for different 
types of waves, including light and sound. Recently, di- 
rect observations of the Anderson localization by disorder 
and of a localization transition by quasiperiodic po- 
tentials 3 of quasi one-dimensional (ID) matter waves of 
ultracold atoms were reported. These experiments pave 
the way to the observation of new phenomena and shed 
new light on long standing problems, amongst which the 
question of possible Anderson localization in presence of 
interactions. 

In the present paper we consider the case of an atomic 
vapor described as a weakly interacting Bose gas in the 
presence of a weak disorder (what is meant by "weak" 
here will be made quantitative in Sec. [TTJ) . In this con- 
figuration it has been shown theoretically in Refs. 0, Q 
and supported by numerical simulations in Ref. Q that 
a small amount of disorder does not drastically alter 
the equilibrium properties of the system, but merely de- 
creases the condensate and superfiuid fractions. Further- 
more, even in the ID limit considered in the present work, 
it has been experimentally demonstrated in Refs. @, [1| 
that one can observe global phase coherence in the pres- 
ence of disorder and remain far from, say, the Bose glass 
phase originally proposed by Giamarchi and Schulz and 
Fisher et al. [3G3]. 

Here we are interested in transport properties. Specif- 
ically, we study a quasi ID, weakly interacting Bose Ein- 
stein condensate (BEC), propagating through a disor- 
dered potential. In this context, localization has been 
theoretically studied mainly for effective attractive in- 
teractions (see, e.g., [ll| and references therein), with 



less attention on the repulsive case we consider here (see, 
however, Refs. [H, [HI]). In the absence of an external 
potential, (repulsive) interactions make the system su- 
perfiuid and introduce a new characteristic speed in the 
system, the speed of sound c. As mentioned above, when 
the speed V of the BEC relative to the external poten- 
tial tends to zero, the addition of a weak random po- 
tential preserves superfluidity, although with a reduced 
superfiuid fraction. What happens as V increases ? This 
question was investigated in a previous publication [Tij . 
where the disordered potential, of length L, was modeled 
by a series of randomly located delta peaks. For small 
velocities V/c <C 1, perturbation theory shows that the 
superfluidity is preserved, e.g., the flow is dissipationless 
and with a perfect transmission. In contrast, in the high 
speed limit 7/c> 1, where the kinetic energy dominates 
over the interaction energy, the transport properties of 
the BEC are deeply altered, and tend to those of the 
non-interacting gas, displaying an exponential damping 
of the transmission with length L, a behavior character- 
istic of the strong Anderson localization. Thus, two lim- 
iting cases of stationary flow have been identified [ijj], 
with contrasting transport properties: superfluidity in 
the deep subsonic regime, and Anderson localization in 
the deep supersonic one. In between, in the region V ~ c 
where both interaction and kinetic energies are impor- 
tant, it was shown that stationary scattering solutions 
do not exist: one reaches a regime of time-dependent 
flows with more or less (depending on the speed) com- 
plex density excitations. The range of speeds around 
c where this phenomenon is observed increases as the 
length L increases. The different types of existing flows 
are summarized in Fig. [TJ 

In the present study we concentrate on the supersonic 
stationary region of the phase diagram [gray (light blue 
online) V/c > 1 region in Fig.[Tj. In this domain we pro- 
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FIG. 1: (Color online) Transport of a quasi ID BEC with 
velocity V through a disordered potential Us consisting in a 
series of uncorrelated delta peaks extending over a domain of 
size L [cf. Eq. (|2ip and the discussion in Sec. IV C] , Dark 
region: time dependent flow; light gray (light blue online) 
regions: stationary flow. In the supersonic case, the yellow 
solid line corresponds to the threshold L* between these two 
domains as determined from Eq. (|81[1 . The blue dashed line 
is the localization length Li oc (|87|) . The supersonic region 
below L\ oc denoted as "ballistic" corresponds to the region 
where the perturbation theory of Sec. IV Al applies. Note the 
enlarged scale for V/c G [0, 1]. 



vide analytical and numerical evidence of Anderson local- 
ization in the presence of interaction for different types 
of disorder. We compute analytically the interaction- 
dependent localization length as well as the correspond- 
ing distribution of transmission coefficients. We also ex- 
plain the disappearance of the supersonic stationary flow 
observed at a given velocity for increasing length of the 
disordered sample. This onset of time dependence is an 
important qualitative effect revealed by our study. We 
show that it is directly connected to interaction effects 
and provide an analytical estimate of the length L* of the 
disordered region above which most of the realizations of 
the random potential lead to time dependent flows (see 
Fig. [J). 

The paper is organized as follows. In Sec. |TT]we present 
the model and identify its range of validity. In Sec. IIIII 
we take some time to properly define the transmission co- 
efficient of a Bose-condensed beam over an obstacle. In 
HVl we introduce the different types of disordered poten- 
tials studied in the present work. In Sec. [V] we present 
analytical and numerical results showing that Anderson 
localization is indeed possible in the supersonic regime. 
We consider the three possible supersonic regimes: per- 
turbative (Sec. IV A[) . Anderson localized (Sec. IV B[) and 
onset of time dependence (Sec. IV C[) . In Sec. |VT]we dis- 
cuss experimental strategies and possible signatures for 
the observation of Anderson localization in an interacting 
Bose-Einstein condensate. Finally we present our conclu- 
sions in Sec. IVIII Some technical points are given in the 
appendixes. In Appendix |A1 we derive the probability dis- 



tribution of transmission in a special case (perturbative 
regime and correlated Gaussian potential) . In Appendix 
[B]we present the derivation of the Fokker-Plank equation 
([55)1 for the distribution of the transmission coefficients. 

II. MODEL 

We study here the transport properties of a quasi- 
one-dimensional (ID) Bose-Einstein condensate formed 
of particles of mass to, experiencing a repulsive effective 
interaction (characterized by the 3D s-wave scattering 
length a > 0) , in the presence of an obstacle represented 
by the external potential U. The potential is not nec- 
essarily disordered at this point, the only restriction we 
impose throughout the present work is that it should have 
a finite extent, i.e., U{x) — > when x — ► ±oo. The con- 
figuration we consider corresponds to the "ID mean field 
regime" [l5[ (see also the discussion in Ref. [Hj]), where 
the system is described by a ID order parameter ij](x,t) 
depending on a single spatial variable: the coordinate 
x along the direction of propagation. if>{x, t) obeys the 
nonlinear Schrodinger equation 

In all the present work we choose to work in the "labora- 
tory frame" where the condensate is initially at rest. Eq. 
([T]) describes its ID dynamics in the presence of an ob- 
stacle moving at constant velocity V in this frame, which 
corresponds to the experimental situation where an ob- 
stacle is swept through a condensate initially at rest, see 
e.g. Refs. [13, [H, EH ■ On the theoretical side, one should 
imagine that, from an initial static configuration where 
the condensate is at rest with U = 0, the potential in- 
tensity and speed have been slowly ramped up to a point 
where the condensate dynamics is described by Eq. |T]). 
We choose V > 0, this corresponds to a potential moving 
from left to right in the laboratory frame. 

The reduction of the motion of the condensate to a 
single spatial dimension is typically achieved through a 
transverse harmonic confining potential of pulsation ujj_ . 
We choose a normalization such that n(x,t) = \i/j(x, t)\ 2 
is the linear density of the condensate. In this case, the 
interaction amongst particles results in E q. (TO in the 
nonlinear term gl^l 2 , with g = 2 hui± a [2(il l2lLl22l|. 

In the stationary regime, where the flow is time- 
independent in the frame moving with the potential, ip 
depends on x and t only through the variable X = x — Vt. 
The appropriate boundary condition is ip(X — > — oo) = 
y/nQ (where no is a constant) (see [22j | and the discussion 
in Section fill Al below) . The condensate is then charac- 
terized by a chemical potential fi = g no, a speed of sound 
c = (gno/m) 1 / 2 and a healing length £ = H/(mc). 

It is customary to characterize the transverse con- 
finement via the "harmonic oscillator length" aj_ = 
(H/muj^) 1 / 2 . With n\ denoting a typical order of mag- 
nitude of n(x,t), the ID mean field regime in which Eq. 
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([T]) is valid corresponds to a density range such that 

{a/a ± ) 2 < rixa< 1 . (2) 

In this domain the wave function of the condensate can 
be factorized in a transverse and a longitudinal part 
[20L I2H [13] • The transverse wave function is Gaus- 
sian (this is ensured by the condition m a -C 1), the 
longitudinal one is of the form ij)(x,t) exp{— i/xt/ft} and 
ip(x, t) satisfies Eq. (JJ) [U[I3]. The left-hand side (l.h.s.) 
inequality in Eq. ([2]) prevents the system to enter in 
the Tonks- Girardeau regime. More precisely, a general 
analysis of ID Bose gas shows that at zero tempera- 
ture no BEC is possible [23| . This results in a alge- 
braic decrease of the one body density matrix monitored 
by phase fluctuations occurring over a phase-coherence 
length = £exp{7r a ± (n 1 /2a) 1 ^ 2 } Hence the 

results obtained using Eq. fl} are valid if they describe 
structures with a characteristic length scale smaller than 
Lff,. The l.h.s. inequality in Eq. ^ ensures that 
is exponentially large compared to the healing length. 
If one considers, for instance, 87 Rb or 23 Na atoms in 
a guide with a transverse confinement characterized by 
u>± = 2tt x 500 Hz, the ratio a/a± is roughly of order 
10~ 2 and restriction ([2]) still allows the density to vary 
over four orders of magnitude. 

Even if the mean field approach is legitimate in ID, 
the effects of disorder have to be taken into account with 
some care. It may well be that the introduction of a 
disordered potential U in Eq. (jTJ) modifies the properties 
of the ground state. This is indeed the case as shown 
in Refs. 0, Q : a disordered potential decreases the 
condensate and the superfluid fraction, but the effects 
are weak provided the intensity of the disorder remains 
weak (see Ref. 26] for an extension to finite tempera- 
ture) . More precisely, in the case of a disorder formed by 
randomly spaced delta impurities with density n$ (see 
Sec. lIV Al) one can show [lj, [27| that, in the dilute impu- 
rity limit, at V = the non-superfluid fraction (normal 
part) is proportional to ns£,{£,/b) 2 [the notations are those 
of Eq. (|2"Tj) ] and thus remains small provided the dimen- 
sionless coefficient (£/&) is small (weak disorder limit). 
At finite V , the normal fraction is multiplied by a fac- 
tor [1 — (V7c) 2 ]~ 3 / 2 (see Ref. (HI), which diverges when 
V = c. One thus expects the mean field approach to 
fail near the region V ~ c of Fig. [TJ This is supported 
by the numerical results presented in [28[. Hence, in the 
center of the time-dependent region of Fig. [1] we can- 
not trust the results obtained from Eq. ([T]). However, 
far from this region, the ID mean field approach is ex- 
pected to be valid even in presence of (weak) disorder, as 
experimentally demonstrated in Refs. @, Q • 

III. DEFINITION OF THE TRANSMISSION 

In the present work we characterize the localization 
properties of the condensate in the random potential by 
studying the transmission coefficient. Eq. {!]) being non 



linear, the definition of transmission and reflection coef- 
ficients needs to be treated with special care. This is the 
purpose of the present section where we first define the 
stationary regime (Section IIII A[) and then the transmis- 
sion coefficient within this regime (Section IIII B[) . 



A. Stationary regime 

It is customary to perform a Madelung transforma- 
tion and to write ip(x,t) = \J n(x, i) exp{i S(x, t)} where 
n(x,t) is the density and hd x S/m — v(x,t) the local 
velocity. From fl]) one can check that they verify the 
continuity equation 

d t n + d x (nv) = . (3) 

The stationary regime is defined as the regime where the 
system is at rest in the frame moving with the obstacle. 
In this case, in the laboratory frame ip, S, n and v are 
time dependent, but they depend on x and t only through 
the variable X = x — Vt. It is then possible to get a first 
integral of ^ under the form 

^(si-')-""' < 4 » 

In the case of subsonic (V < c) and stationary mo- 
tion, the flow is superfluid and the order parameter is 
only affected in the vicinity of the obstacle, with n(X — > 
±00) = n and v(X -> ±00) = [HHH]. 

For V > c, a regime of stationary flow also exists but 
in this case the obstacle induces density oscillations with 
a pattern stationary in its rest frame [22| . This means 
that in the laboratory frame the phase velocity of these 
waves is identical to the velocity V of the obstacle. On 
the other hand, the energy transferred from the obstacle 
to the fluid propagates with the group velocity, which 
in the case of Bogoliubov excitations is greater than the 
phase velocity, i.e. - as just argued - than V. As a 
consequence, radiation conditions require that the wake 
is always located ahead of the obstacle, i.e., upstream, 
with no lon g-ra nge perturbation of the fluid on the down- 
stream side [22|, |3(| • This means that in this case the flow 
far in the downstream region remains unperturbed, with 
n(X — ► — 00) = tiq and v(X — > —00) = 0. The two pos- 
sible stationary configurations (subsonic and supersonic) 
are represented in Fig. [21 

Hence, in any stationary configuration (subsonic or su- 
personic), the above reasoning fixes the integration con- 
stant in the right hand side (r.h.s.) of Eq. ^ to its value 
at X — > — 00, i.e., — uqV. 

In the stationary regime one gets from Eqs. (JTJ) and 
© 
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FIG. 2: Schematic representation of the typical density pro- 
files. The upper plot corresponds to a subsonic stationary 
profile, while the lower one corresponds to a supersonic sta- 
tionary profile. The potential moves from left to right, and 
the upstream (downstream) region thus corresponds to the 
region X — * +oo (X — > — oo). In both plots the potential is 
represented by a thick solid line (hatched down to zero) and 
the density profile is represented by a thin solid line. 



where A(X) = y/n(X)/no and 



W(A) = ^(A>-1) 



V 2 

V z - c z A z - — 
A 2 



(6) 



B. Transmission coefficient 

In this section we restrict the analysis to the station- 
ary regime of section MI Al and define the transmis- 
sion of the condensate through the obstacle represented 
by a potential U (not necessarily disordered) verifying 
U(\x\ -> oo) = 0. 

As the wave equation ([1]) is nonlinear, one cannot, in 
general, properly define reflection and transmission coef- 
ficients, since it is generally not possible to disentangle 
incoming and reflected waves in the nonlinear flow up- 
stream the obstacle. However, following a procedure de- 
vised in Ref. [3l[ (see also [32]), we will show that one 
can define a transmission and a reflection coefficient in 
the limit of small nonlinearity as well as in the limit of 
weak reflection and arbitrary nonlinearity. 

Outside the scattering region, U(X) = and one can 
get a first integral of Eq. |5j) under the form 



h 2 

2m 



dA\ 



— \ax) +W(A)=E*, when X - ±oo , (7) 



which defines the "free" asymptotic density profiles. E^ 
in Eq. ((7]) are integration constants. The boundary con- 
dition discussed in the previous section imposes A = 1 
and dA/dX = when X — > — oo. This fixes the value 
E~, = 0. The value of Et at +oo has to be determined 



FIG. 3: (Color online) W as a function of A = \ijj\/y/no 
(drawn for V/c = 4). The fictitious particle has a classical en- 
ergy when X — ► +oo. The (red) dashed line corresponds 
to an approximation of 14^(^4) by h 2 K 2 (A — l/A) 2 /(2m), ob- 
tained by keeping only the first term in expansion (|13[) . 



by the integration of Eq. fl]) (cf. Ref. [22]). Eq. ex- 
presses the energy conservation for a fictitious classical 
particle with "mass" ti 2 /m, "position" A and "time" X, 
evolving in a potential W (whose typical shape is dis- 
played in Fig. [3]) • This type of analysis is common in the 
study of nonlinear equations such as Eq. ([T|), see e.g., 
the review [33[ (the first time we found it used is in Ref. 
[l^]). It is employed here as a convenient tool for get- 
ting intuition about the behavior of the solution of the 
Gross-Pitaevskii equation (see below). 

From now on we restrict to the supersonic stationary 
regime where an imperfect transmission occurs (in the 
subsonic stationary regime one has perfect transmission) . 
In this case the fictitious particle is initially (i.e., when 
X — > — oo) at rest at the bottom of the potential W 



with E, 



0. 



on the value of E^ 



The behavior of A for X — ► +oo depends 
A stationary solution exists only 
if A(X —>■ +oo) remains bounded, i.e., if E^ < W(Ai) 
(where A\ corresponds to the local maximum of W , see 
Fig. [3]). In this case, the asymptotic behavior of A(X) 
corresponds to oscillations between the values ^4 m i n and 
A max defined in Fig. 3. 

For future references we note that 14^(^4) is zero when 
A = Aq = 1 and when A = V/c, and that the derivative 
diy/dA is zero when A = Aq = 1 and when A = A\, 
with 



(8) 




At large velocity, when V ^> c, one has 
V „ / 1 



and 



Ai = 



W(A X ) 



V / 2( 



O 



V/c 



8 c 2 



:v 2 



(9) 



(10) 
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Writing A 2 (X) = pjX) = 1 + 5p(X) we now argue 
that, following Ref. [3l| . one can write a perturbative 
version of Q in a limit where 



\Sp(X)\ « 



V 2 



(11) 



We emphasize that restriction (fTTj) corresponds to \Sp\ -C 
1 (i.e., small oscillations) only when V is of order of c or 
smaller. The approach developed below is however able 
to tackle large relative density oscillations (\6p\ ^> 1) at 
large velocities (F>c) (35|. In this sense it will allow 
us to penetrate in the non perturbative regime where the 
upstream density oscillations are large and the transmis- 
sion is low. 

Using the variable p, we write Eq. 10 in the upstream 
region (X — > +oo) as 



2m 



dp_ 

dX 



8F(p) = 8pE i 



cl ' 



(12) 



where F(p) = pW(A = -y/p). A simple limited expansion 
around p = 1 yields 



Hp) 



2„2 



H K 



2m v 2 
where 8p(X) = p(X) — 1 and 



« = 2|v a -.?| 1/a 
h 1 1 



and 



(13) 



K £ . (14) 



The second term in the r.h.s. of expansion (fT3")) is small 
compared to the first one precisely in the limit pip . In 
the following we restrict to this regime and neglect the 
second term of the r.h.s. of p3[) . This corresponds to ap- 
proximating the exact W (A) by the (red online) dashed 
line in Fig. [3] and to write Eq. (|12p under the form 



+ An 2 (8p) 2 = 16k 2 A(1 + Sp) 



where the dimcnsionlcss parameter A is defined by 

mE+ 

2h 2 K 2 ' 



(15) 



(16) 



The solution of Eq. (TJ5]) is 

n(Jf) 



no 



p(X) = 1 + 2A + 2Acos(2kX + 6») , (17) 



where 



A= VA2 + A 



(18) 



and # is a phase factor. We recall that Eq. (|17|) describes 
the density oscillations in the upstream region X — ► +oo. 
These oscillations can be described as the sum of incident 
and reflected waves (V'inc and ^ re f) of the form 



V'inc(A) = ^0(1 + A) exp(-mX) 
i/j ie i(X) = \J noX exp(incX + i0) . 



(19) 



This analysis allows one to determine the reflection and 
the transmission coefficients as 



R = 



\Aef\ 2 
IVincP 



A 



1 + A 



T=l-R= 



1 



1 + A 



(20) 



Of course the sum of the incident tpi nc and the reflected 
ij) Te f waves (|19p is an approximate solution of the nonlin- 
ear Schodinger equation |T]) which is only valid in regime 
(jlip . i.e., in the regime of arbitrary interaction and small 
transmission (A <C 1), or in the regime of arbitrary trans- 
mission and small interaction (F»c). 



IV. DIFFERENT TYPES OF DISORDER 

Up to this point we presented a theory valid for any 
potential of finite extent. From now on we concentrate on 
the particular case of random potentials. We denote U(x) 
an arbitrary random potential, and use a subscript when 
dealing with one of the particular cases defined below. 



A. Potential formed by a series of 5 peaks 

The first potential of interest, analyzed in Ref. [14| , is 
a series of N randomly located identical delta peaks of 
the form : 



U 5 (x) 



(21) 



The intensity of the peaks is measured by the (non ran- 
dom) positive quantity b. The scatterers have random 
uncorrelated positions = X\ < X2 < x%..., with mean 
density ng and average separation lg — 1 /rig. Hence the 
potential extends over a mean length L = (N — l)lg. 

Denoting henceforth the disorder average by (..), for 
x and x' inside the disordered region one gets the mean 
value 



{Us(x)) = 



h 2 ng 
m b 



(22) 



and the irreducible two-point correlation function 

(Ug{x)Ug(x')) - (Ug) 2 = a (H 2 /m) 2 6(x - x') , (23) 
where 

(24) 



b 2 ' 



B. Correlated Gaussian potential 

Another commonly used model of disorder is provided 
by Gaussian random processes with zero average. We 
consider here potentials which are non zero only over a 
region of finite extent (with typical size L) and generate 
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them in the following way (see, e.g., Chap. 5 of Ref. [36[ 
and references therein): let's consider a Gaussian white 
noise i](x) extending over all the real axis, with zero mean 
and second moment (r](x)r](y)) = 5(x — y). Then for a 
given function w(x) one defines 



U g (x) 



w{x - y) rj(y) dy , 



(25) 



where a is a parameter characterizing the disorder and 
whose meaning is explained below. If w were a delta func- 
tion, then U g would be a Gaussian white noise over [0, L] 
(and zero everywhere else). The actual function w(x) has 
a finite extension, and this induces finite correlations in 
the disordered potential. 

From l]25p it is clear that (U g ) = 0. If the domain of 
integration in the r.h.s. of (|25p were extended to all M., 



(26) 



U g would have a Gaussian distribution 



exp[ 



25? 



V2ttS 2 



where 



Y? = a (h 2 /m f [ w 2 (x)dx . (27) 



Defining the correlation function C as 

(Ug(x)Ug(x')) - (Ug) 2 = V (ft 2 / ' C (x ~ X ' ) 

one would get in this case 



C(x) = dyw(x + y) w(y) , 
with a Fourier transform 

C(q) = [dxC(x)e- i ** = \w(q)\ 2 , 
it 

where w is the Fourier transform of w. 

Imposing here the normalization condition 

w(x) dx = 1 , 



(28) 



(29) 



(30) 



(31) 



leads to a two-point correlation function (|28|) which is - 
as in Eq. (f2"3")l - of the form of er(ft 2 /m) 2 multiplied by a 



function whose integral over x equals unity \C(x) in Q28[) 
instead of 5(x) in |23p]. Thus, with definition (|2"5")1 and 
normalization (|3lTl . cr plays for disorder ([25]) the same 
role as ng/b 2 [Eq. (|24p] for disorder (|2ip : it characterizes 
the amplitude of the fluctuations of the potential. The 
typical extent of w(x) will in turn characterize the range 
of the correlations. 

Since U g as given by Eq. (|25p is typically non zero only 
over a region of finite extent, Eqs. ([28]) and ([29]) are only 
correct if x and x' are inside this region. More precisely, 
they should be in this region, at a distance from or 
L larger than the typical extent £ c of the function w. 



In the following we always consider the case where L is 
very large compared to £ c (otherwise one could simply 
not speak of a disordered region) and it is clear that the 
characteristics of the disorder are properly defined only 
inside the disordered region. 

We consider two special cases of correlation corre- 
sponding to different forms of w: a Lorentzian 



, x 1 4/2 
w L {x) = - 



it (4/2) 2 + x 2 



and a Gaussian 

w G {x) 



1 l X 

exp 



(32) 



(33) 



We denote the corresponding potentials by U L and U G 
For the correlation functions one gets respectively 



C L {x) 



4/tt 
£ 2 + x 2 



C L (q)=e~ 



:\Q\ 



(34) 



and 



C G (x) 



-x 2 /(2£ 2 c ) 



V^e 2 



, c G ( q )=e-" 2 ^ 2 . (35) 



In both cases l c is the typical correlation radius. 

The choice of a Lorentzian correlated disordered poten- 
tial originates from experimental and theoretical results 
in the case of micro- fabricated guides. In this type of 
setting, the atoms are magnetically guided over a chip 
[37l | . Unavoidable imperfections and irregularities in the 
design of the circuit induce fluctuations in the current 
which, in turn, result in a random contribution to the 
magnetic field used for guiding the atoms. Thus the 
potential seen by the atoms has a random component 
which is typically Lorentzian correlated, with a correla- 
tion length £ c which decreases when the distance between 
the guide and the chip increases [HI, HE HH, Ell ■ The 
Gaussian correlated potential U G is more academic but, 
by comparison with the results obtained with U L it al- 
lows one to check what is really specific to the Lorentzian 
case, and what is a mere effect of finite correlation length. 



C. Speckle potential 

Another experimentally relevant type of disorder is 
the so called speckle potential which is generated by an 
optical speckle field produced b y a laser beam passing 
through a diffusing plate [HHH,!!!]- The corresponding 
potential will be denoted by U s and may be mathemati- 
cally generated as follows [44j |: 



U s (x) 



L 

w s (x 



y) [vi(y) + iv2(y)]dy 



(36) 



where 771 and 772 are two independent Gaussian white 
noise processes of zero mean with (r) a (x)r) a / (x')) = 8(x — 
x')5 aa i (a and a 1 = 1 or 2). 
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Here also we characterize the disorder by studying its 
statistical properties in the limit where the domain of 
integration in the r.h.s. of |36|) is extended to all R. In 
this case one gets 



nu a ) 



exp(-^r) 
2S 2 



(37) 



where S is given by formula (|27p (replacing w by w s ). 
This yields (U s ) — 2Y? and the correlation function de- 
fined in Eq. (T3S|) reads here 



C s (x-x') 



1 



a (h 2 /m) 2 



[{U s (x)U s (x'))-{U s ) 2 ] 

2 (38) 



dyw s (x - x' + y)w s (y) 



=4 



Contrarily to the choice (f31j) . w s should not be nor- 
malized to unity here because - from Eq. ([36]) - this is 
homogeneity-wise impossible. Instead, the choice 



1/4 sin (^) 
4 IT 6 J X 



(39) 



corresponds to the typical experimental situations [42 
and leads to a correlation function 



C 3 (X) 



(40) 



whose integral over x equals unity and whose Fourier 
transform is 

1-M4/2 if \q\t e <2, 
Cs{q) = { (41) 
otherwise . 



Hence, definition (|36| and choice (|39|) correspond here 
also to characterizing the amplitude of the disorder's fluc- 
tuations by the parameter a and the range of the corre- 
lations by £ c . 



V. SUPERSONIC STATIONARY REGIME 

As explained in Ref. [l4| , and recalled in the Introduc- 
tion, Anderson localization in a weakly repulsive Bosc- 
Einstein condensate is only possible in the supersonic 
regime (cf. Fig. [1} which we consider now. In the 
present section we first analyze the transmission across a 
short disordered sample, in which case perturbation the- 
ory is applicable (Sec. IV A[) . We then turn to generic 
non-perturbative configurations (Sec. IV Bp where An- 
derson localization is expected. In this regime we obtain 
evidences of the occurrence of Anderson localization in 
the presence of interaction. Finally we discuss the up- 
per limit of the localized regime and the onset of time 
dependent flows for long disordered samples (Sec. IV C|) . 



A. Perturbation theory (A <C 1) 

In the supersonic stationary regime, simple perturba- 
tion theory yields n(x, t) — no + 5n(X) where [22J 



Sn(X) 



2 m no 
h 2 k 



x 



dyU(y) sm[2 K (X-y)} 



(42) 



and k is given by Eq. (|14[) . Perturbation theory always 
predicts a stationary density profile. This is certainly 
wrong when V is close to c (cf. Fig. []}, but in this 
case k gets very small and one precisely goes out of the 
domain of validity of perturbation theory (Sn as given by 
(|42| is no longer small compared to no). 

Far ahead of the obstacle (in a region where X — L is 
larger than £ c and (|4"2")l gives 



5n(X) _ 2m 



"o 



(43) 



The perturbative regime in which Eqs. (|42|43|) are valid 
is also the one where the constant A in (|16I17[) is small 
compared to unity. This can be inferred from the com- 
parison of ([4"3"|) and (TIT)) which indeed shows that A C 1 
in the regime where (|43|) holds and that, in this case, 
\/A ~ m\U(2K)\/(H 2 n). The corresponding reflection co- 
efficient can then be obtained from ([20]) . yielding 



R~ A: 



h 4 K 2 



\U(2n)\ 2 



(44) 



From (|4"4"]) it is clear that the average reflection coefficient 



TD 

W = <A) = ^<|C/(2k)| 2 )<<1. (45) 

Furthermore, one can show that the corresponding prob- 
ability distribution of the reflection coefficient is Poisso- 
nian with 



P(R) = -±rexp(-R/(R)). 



(46) 



Note that for properly normalizing this probability dis- 
tribution for R S [0,1], one should include in the pref- 
actor of the r.h.s. of (|46|) a correcting term of order 
cxp(— l/(R)) which can be safely neglected in the limit 

6S}. 

We give in Appendix [5] a demonstration of result ([46]) 
for the special case of a correlated Gaussian potential U g 
of type (|25|) . Below, we show that the same result holds 
for a potential Us of type (|2"Tj) [see Eq. (fB"5|) in section 
IV B| , and we checked numerically that it is also the case 
for the speckle potential U s (|3"6"j) (cf Fig. 2]). In all these 
cases, the average reflection coefficient reads (up to the 
above discussed exponentially small correction) 



(R)=L/L 1oc {k) , 



(47) 




FIG. 4: Probability distribution P(T) for the transmission 
coefficient T in a potential Us with a = 3.14 u 2 £, £ c = 0.1£ 
and L = 50 £ moving at velocity V = 7c in a condensate of 
initial constant density no £ = 1. The corresponding localiza- 
tion length is Li oc = 500 £. The histogram corresponds to a 
statistical analysis of the results of the numerical solution of 
Eq. |T| for 10000 different random potentials. The solid line 
is the perturbative result (|49[) . 



where 



^Ioc(k) 



n 2 /a 
C(2k) 



(48) 



We recall that the function C depends on the type of 
disorder considered. For a potential of form (|21[) one has 
C$ = 1, for the other potentials considered in this work 
it is given bv p4|) . ([35]) and (j41"]l . 

Concomitantly to distribution (|4"rj|) of reflection coeffi- 
cients one gets for the transmission 



P( T ) = ^exp(-(l-T)^ 

Li 1 Lj 



From (gSJ [or one gets 



(T) 



^Ioc(k) 



(49) 



(50) 



The perturbative approach holds when (R) <C 1 i.e., 
when L <C £i oc - This corresponds to the region which 
is denoted as "ballistic" in Fig. Q] (45|. Its accuracy is 
shown for L/L\ oc = 0.1 in Fig. Q] for a speckle potential 
U s of type ([55)1 (we also checked this prediction for the 
potentials Us and U G , with also excellent results). 

At this stage, L\ oc is simply a notation for expression 
(l4"8"| . but it will be shown to be the actual localization 
length of the matter wave in a disordered potential (in 
section |VB|) . 

The results derived here also hold for a noninteracting 
gas, obtained by taking the limit g — > in Eq. {!]), in 



FIG. 5: (Color online): Average transmission as a function of 
L for a potential Us (characterized by ns£, = 0.5 and £/& = 
0.1). In the interacting case V = 3c and L\ oc (k) = 1600 £. 
The non-interacting case is drawn for the same velocity and 
corresponds to a value L\ oc (k) = |Li oc (k) = 1800 £. In both 
cases the dashed line is the analytical result (|50p and the 
solid line corresponds to a statistical analysis of the results of 
the numerical solution of Eq. ^ for 15000 different random 
potentials. The departure of the numerical results from the 
dashed lines occurs when the systems leaves the perturbative 
regime. 



which case c = and k — k. Eq. (|48p. with k replaced by 
k, then coincides with the Antsygina-Pastur-Slyusarev 
formula for the localization length [46l . [47J and the distri- 
bution of transmissions (|49[) holds, with L\ oc = L\ oc {k). 

In the present work, those formulas are modified to 
include interactions. The generalization simply con- 
sists in replacing the wave vector k = mV/h by n — 



m(V 2 



') ' /%. This replacement has, as an important 



physical consequence, the effect of diminishing, at a given 
speed V, the localization length (there is an effective re- 
duction of the available kinetic energy by the repulsive 
interactions). For instance, in the case of a potential Us, 
since L\ oc (k) oc k 2 , there is a relative difference c 2 /V 2 
between Li oc (k) and Li oc (k), that is 11 % for V = 3c. 
This is illustrated in Fig. El which displays the average 
(T) as a function of L for a disorder Us of type (p?Tj) , with 
and without interactions. 



B. Non perturbative approach 

When the size L of the sample is large compared to 
the value L\ oc determined in Section TV Al the perturba- 
tive approach fails. We now propose a non perturbative 
method allowing to treat both the regimes L < L\ oc and 
L > L\ oc and showing that L\ oc , as defined in Eq. (|48]) . 
is indeed the localization length in the presence of inter- 
actions. 

Within the framework of the non perturbative ap- 
proach, we are able to provide approximate analytical 
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E 



(i) 



cl 



E 



(2) 



x/% 



FIG. 6: (Color online) Upper panel: W as a function of A 
(drawn for V/c — 4). For X < xi = 0, the fictitious particle 
is initially at rest at the bottom of potential W with E t 



(0) _ 

cl — 
(n-l) 
cl 



0. The value of the classical energy changes from E t 

(n) 

to E^i at each impurity x n . The lower panel displays the 
corresponding oscillations of A(X), with two impurities at 
xi = and £2 = 4.7 £ (their position is indicated by vertical 
dashed lines). 



results in the case of the model disorder potential (|2ip . 
This potential being zero between two impurities, one can 
write a series of first integrals of Eq. ([5]) in each segment 
]x n ,x n +i[ as follows : 



tf_ fdAV 
2m YdXJ 



W[A(X)] = E, 



(n) 



(51) 



In the region X < x\ = the integration constant E~ x of 
Eq. f7J is denoted as E^ in (j5"Tj) (taking x = —00) and 
is zero, whereas in the region X > xn one has E^ = E^ 
(and .tjv+i = +00). 

From ||SJ) it is a simple matter to show that the match- 
ing condition of the density at impurity position x n is 



A'(x+)-A'(x-) = -A{x n ), 



(52) 



where A'(x~) L4'(x + )] denotes the limit of the derivative 
dA/dX at the left [at the right] of x n . Relation ([52]) 
between the derivatives of the amplitude results [from 
Eq. ([ST) ] in a relation between the classical energies: 



E. 



(n) 



E 



(n-l) 



2h 2 

mb 2 



A(x n )[bA'(x-) + A(x n )] . (53) 



Hence, Eq. ([ST]) allows to draw a classical analogous 
of the solution of the nonlinear Schodinger equation in 
the presence of potential (|2"Tj) formed by a series of delta 
peaks: the fictitious classical particle defined in Sec. 
IIIIBI evolves in the potential W and experiences kicks 
at "times" x n . Each kick changes the "energy" accord- 
ing to (fS"3"]) . as illustrated in Fig. [Sj The key point in the 



remaining of this section will be to derive the probabil- 
ity distribution of E^ = E^ which then directly allows 
one to get the distribution of A's and of the transmission 
coefficients [through Eqs. ((TBJ) and ([50])]. 
Let us introduce the quantities 



" ~ 2K 2 k 2 



and A„ = yj (A„) 2 + A„ 



(54) 



The parameters A and A defined in (|16p and (fT5|) are 
related to the ones of Eq. (f5"4"]) by A = Xn and A = Ajy 
(i.e., A is the last of A„'s; the same holds for A). Denoting 
by On-i the value of the phase [appearing m p"?]) ] for 
X £]x n —i, x n [, one gets in this domain [the derivation is 
exactly the same as for Eq. fP7|) ] 

A 2 (X) = 1 + 2A„_! + 2A„_x cos (2kX + 6>„_i) , (55) 



and one can rewrite Eq. (|53|) as 

2h 2 



E 



{n) 



E 



(n-l) 



lb 2 



[l + 2A„_i] 



Ah 2 

-^An-lV K 2 b 2 + 1 Cn-1 , 

mb z 



where 



Cn-i = cos [2 KX„ 



+ tan^ 1 (nb)] 



(56) 



(57) 



Using definition (|54|1 one can rewrite Eq. (|56p in terms 
of the parameter A n as 



A n — A n _ 



1 + 2 A„_i 2 A„_i 



K 2 b 2 



K 2 b 2 



V* 2 b 2 + 1 Cn-i • (58) 



Eqs. (|56I58|) are valid provided (fTT|) holds, i.e., pro- 
vided E { ™~ 1] < W(At), which reads 



E 



(n-l) 



h 2 n 2 



< 



(k 2 C ), ».e., A„_i« K ^ ■ (59) 



In the following we also impose the condition 
k6> 1 . 



(60) 



Precisely, we neglect all the quantities of order 1/(k 3 & 3 ). 
This is an important technical point. It facilitates the 
analysis by allowing one to get simple formulas as we 
now illustrate in the perturbative case: Eq. (|55|) allows 
for instance to compute the average value of the reflection 
coefficient in the perturbative regime (as already done in 
Section IV A[ Eq. (|47)l ). In this regime, additionally to 
condition (|59|) one has A„ C 1. Then Eq. (|58|) implies at 
leading order (A„) = (An-i) + 1/(k 2 6 2 ) which, together 
with the initial condition Aq = 0, leads immediately to 



(A) 



K 2 b 2 



(N) 



1 L 



n 2 b 2 l s ' 



(61) 



which is identical to result (|47p in the case of a potential 
Us for which Li oc (k) — K 2 /a — n 2 b 2 Is [cf. Eq. 
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Let us now proceed and consider the generic non per- 
turbative regime where A„ may become large compared 
to unity and where (jTTJ) and (|59|) are still valid. Taking 
the condition (|60p into account, Eq. ([55)) reads 



A n — A n - 



l+2A ra - 
K 2 b 2 



2A n _ 
k 6 



Cr, 



(62) 



It is natural to assume that the phase of the cosine in the 
r.h.s. of (|57[) is uniformly distributed in [— 7r, 7t] and inde- 
pendent of the phase at step n — 1. This could be called 
a "phase randomization" approximation. This relies on 
hypothesis (|60j) and on the assumption that there is a 
large number of density oscillations over the (random) 
length between x„_i and x n , i.e., 



(63) 



k{x„ - X n --i) = Klg > 1 



Then, the argument of the cosine in definition (|57|) is 
uniformly distributed, Cn 's are uncorrelated random vari- 
ables, with all the same law characterized by its average 
(Cn) = and variance 



1 



Note that the regimes (j60|) and (l63|) imply that 



ft 2 K 2 h 2 n s 

> - = (Us) , 



2 rn 



m b 



(64) 



(65) 



which in turn implies that the kinetic energy ^mV 2 is 
much larger than (Us); i.e., one is exactly in the An- 
derson regime where the incident kinetic energy is much 
larger than the typical value of the (disordered) poten- 
tial representing the obstacle. Hence a classical particle 
would flow almost unperturbed over the potential but, as 
we shall see, a quantum particle experiences an exponen- 
tially small transmission. 

Let P(X,n)dX be the probability that A„ lies in the 
interval A, A + dA. Going to the continuous limit and 
defining the continuous variable t = n/(K 2 b 2 ) = X/L\ oc 
[where L\ oc = n 2 /a is the parameter (|48| in the case of 
a potential Us] it is shown in Appendix [B] that P(A, t) 
verifies the following Fokker-Planck equation 



dP 
~dt 



d_ 

dX 



dP 

X ^ + X ^8X 



(66) 



Equation (|66|) follows directly from Eq. ([62]) in the regime 
where conditions (fBT)]) and (I63|) hold. It is precisely the 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equation [HI 
for the transmission in a single disordered channel [with 
T = 1/(1 + A)]. Equation (|66|) is sometimes referred to 
as Mel'nikov's equation (after Ref. [4!|) but has a much 
longer history (see the discussion in Refs. pTl. l50l] ) . 

Since before entering the disordered region the particle 
has a classical energy E^' — corresponding to A = 0, 
Eq. ([6"6"1) has to be solved for the initial condition 



limP(A,t) = S+IX) 
t->o 



(67) 



where 6+ is the one-sided delta function: J °° 5 + (A)dA = 
1. In the limit of small t (i.e., in the perturbative regime 
X <C iioc)j A remains small and one can approximate in 
the r.h.s. of (f6"6"|) the term A(A + 1) by A. It is then simple 
to verify that the solution of this approximate equation 
that satisfies (l67l) is 



P(X,t) = 



exp{— X/t} 



t 



for t < 1 



(68) 



This result for the small t solution of the DMPK equation 
has been already obtained in Ref. [5l| (see also the dis- 
cussion in Ref. 52]). The distribution law is exactly 
equivalent to distribution (|46f of the reflection coefficient 
in the perturbative regime and this proves the validity of 
the Poissonian distribution (49) for a potential Us of type 

HUD H. 

In the general case (i.e., for all t > 0) the solution of 
(l6"6l with the initial condition (f6~7| is (see, e.g., Refs. [St], 
[50( and references therein) 



P(X,t) 



-t/4 



V2tt£ 3 



u e 



-« 2 /(4t) 



y cosh(w) — 1 — 2 A 



du , (69) 



where u\ — cosh~ 1 (l + 2A). 

From distribution (|69|) . a- lengthy computation or al- 
ternatively the direct use of the DMPK equation (l66|) [54] 
yields 



>T)= f 
Jo 



dA In 



1 + A 



P(X,t) = -t . 



(70) 



In the large t limit, distribution (|69|) tends to a log- 
normal distribution, i.e., the distribution of the variable 
InT is Gaussian (see Ref. [54] ]) 



P(lnT,i) 



exp{-(t + lnT) 2 /4t} 



for t > 1 . (71) 



From this distribution one gets the correct average 
(InT) = — t, Eq. (T7U|) . and a standard deviation 
[((InT) 2 ) - (InT) 2 ] 1 / 2 = V2t, which is in agreement with 
the exact result in the limit t ^> 1 (54J. At the extrem- 
ity of a sample of length L one has t = L/L\ oc , and the 
distribution (|71[) is the log-normal distribution of trans- 
mission typical for Anderson localization in the regime 
L » Lioc (see, e.g., Ref. [13, HI])- As a side product 
of this analysis, Eqs. ([70]) and (|7Tj) confirm that L\ oc is 
indeed the localization length as was anticipated in the 
notation. 

We have tested the validity of the DMPK approach 
for a Bose-Einstein beam of interacting particles prop- 
agating in a disordered potential Us of type (f2Tj) . The 
numerical results for the probability distribution P{T) 
are compared on Fig. [7] with the DMPK prediction (|69[) . 
The agreement is seen to be excellent. The distribution 
evolves from the Poissonian result (|4"9]) (for low values of 
L/L\ oc ) towards a distribution peaked at low T- values for 
large L/L\ oc . In this latter case one can check that the 
distribution tends to a log-normal by plotting P(lnT). 



11 




0.2 0.4 ^ 0.6 0.8 1 

T 



FIG. 7: (Color online) Probability distribution of the trans- 
mission through a disordered potential Us of type (|2ip (char- 
acterized by £/b = 0.5 and ns£ = 0.5) plotted for differ- 
ent values of the ratio t — L/L\ oc {k) with V = 30c. The 
black solid lines are the DMPK result ([69} and the colored 
histograms correspond to the numerical simulations (50000 
samples used for each value of t). Cases (a), (b), (c) and (d) 
correspond respectively to t = 0.1, 0.5, 1 and 2. The inset 
displays (InT) as a function of t. The thick solid line is ex- 
tracted from numerical simulations and the thin (red) solid 
line is the DMPK prediction (|70[) . They can be distinguished 
only around t ~ 2 as shown in the blowup of the (gray) shaded 
region for 1.9 < t < 2. 



We have also checked the validity of the present ap- 
proach over a sizable range of lengths of disordered re- 
gion and of intensities of disordered potential by plotting 
in the inset of Fig. [7] the average (InT) as a function of 
L/L\ oc . The agreement of the numerical results with the 
DMPK prediction ([70)) is excellent. Note however the 
beginning of a small departure around L/ L\ oc ~ 2; this 
effect will be studied more thoroughly in Sec. IV CI (cf. 
Fig. ©. 

Finally, we discuss numerical results obtained for the 
disordered potentials introduced in Sees. HVBI and llVCI 
Although we do not have an analytical derivation of the 
DMPK equation for these potentials, the numerical re- 
sults indicate a very good quantitative agreement for a 
disordered potential U G and for a speckle potential U s . 
We display the comparison of the numerical data with 
the DMPK predictions for a speckle potential in Figure[5] 
The same agreement is obtained for a Gaussian potential 
U G . Hence, the behavior analytically predicted for the 
potential Us appears to be of general validity, meaning 
that the above defined regime of "phase randomization" 
can probably be extended to correlated potentials, lead- 
ing to a regime of single parameter scaling. However, 
we have noticed that, although showing an overall good 
agreement with the DMPK prediction, the Lorentzian 
correlated potential U L exhibits some deviations in the 
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FIG. 8: (Color online) Probability distribution of the trans- 
mission through a speckle disordered potential Us (character- 
ized by £ c /£ = 0.05 and a = 3.13 (i 2 £) plotted for different 
values of the ratio t = L/L\ oc (k) with V = 13c. Curves 
(a), (b) and (c) corresponds respectively to t — 0.31, 0.52 
and 0.68. For each curve, the black solid line is the DMPK 
result (I69p and the colored histogram is the result of 10000 
numerical simulations. 

tail of the distribution, the details of which will be stud- 
ied elsewhere. 

C. Threshold for the existence of a stationary flow 

In the previous sections the main effect of interaction 
has been shown to be a renormalization of the localiza- 
tion length L\ oc . Interaction induces a modification of the 
wave vector: expression (|48l) for the localization length 
coincides with the noninteracting one but computed for 
an effective interaction-dependent wave vector k given by 
Eq. ([U]), instead of k = mV/h. The repulsive interac- 
tion diminishes the available kinetic energy and therefore 
reduces the localization length with respect to the non- 
interacting case (since k < k). 

We now discuss another, more spectacular, effect of in- 
teractions on the localization properties of a propagating 
BEC on a disordered potential. 

In the previous sections IV Al and IV Bl we completely 
neglected the presence of an upper limit for the classical 
energy E c \, which is given by the local maximum of the ef- 
fective potential W(A), namely £™ ax = W{A{) (see Fig. 
[3]). Trajectories that pass beyond _E™ ax would become 
unstable and develop singularities with infinitely large 
density at X — > oo. In practice this implies, on the level 
of Eq. |T]), that a stationary flow cannot be maintained 
in this case and that the disorder induces time-dependent 
dynamics of the condensate. 

In the vicinity of E™ ax , the density profile of the con- 
densate in between two adjacent scatterers becomes quite 
different from the cosine shape (fTTj) that was derived 
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for weak nonlinearities and/or low density modulations, 
and resembles more to a periodic train of gray solitons 
[22I ]. In a crude approximation, we neglect this com- 
plication and assume that the spatial evolution of the 
density is still given by Eq. (| 1 T[) for all classical ener- 
gies until E c \ = E^ ax . Trajectories that happen to pass 
beyond E™ ax are considered to be "lost", i.e., they do 
no longer contribute to the probability distribution for 
the transmission. This formally amounts to introduc- 
ing a "sink" in the stochastic equation (|62p , namely at 
A = A max = m E™ ax I (2ft 2 k 2 ) . In the corresponding 
Fokker-Planck equation (|66[) . this sink is appropriately 
modeled by imposing the boundary condition 



P(A max ,i) = 



(72) 



As a consequence of this boundary condition, the in- 
tegrated probability distribution J Q max P(A, i)dA is no 
longer conserved, but decreases with increasing t, i.e., 
increasing length L of the disorder region. 

In the following, we show how this affects the DMPK 
predictions of section IV Bl and how the "survival proba- 
bility", i.e., the fraction of trajectories that remain be- 
low this boundary at given length L, can be analyt- 
ically computed in the limit V S> c. In this limit, 
from Eq. (JIDJ) one gets E™ ax ~ mV 4 /(8c 2 ) and thus 
A max ~ V 2 /(16c 2 ) 3> 1. Modifications of the probability 
density P(X, t) due to the presence of the sink appear 
only when the typical value of A is of order A max which, 
as just remarked, is large compared to unity in the case 
V ^> c. In this case P(X,t) is already negligibly small 
around A ~ 1. We therefore make the approximation 
A + 1 ~ A in the Fokker-Planck equation (j6"B"|) , which 
then reads 



dP 
~dt 



d_ 

dX 



A dX 



(73) 



Using, from now on, the probability distribution 
P(\nT,t) for finding a given value of InT at fixed t = 
L/L\ oc , we obtain in this limit 

d d 2 d 

ol p{z ' t) = d^ p{z ^-di p(z ' th (74) 

where we introduce z = — InT. Clearly, the log-normal 
distribution (|TTj) corresponds to a solution of Eq. (fT4")l in 
the absence of any additional boundaries. 

In the presence of the sink, which is imposed by the 
boundary condition P(z max ,t) = with 



ln(A n 



In A n 



In 



V 2 
16c 2 



(75) 



we can straightforwardly find the solution of Eq. (|74|) 
by subtracting from the log-normal distribution (|71[) a 
"mirror" distribution centered at some z > z max (namely 



tribution 

P(z,t) = 



exp 
exp 



t 2z max ) 

it 



(76) 



which is defined for z < z max . Clearly, this distribu- 
tion satisfies the evolution equation (f?4")) as well as the 
boundary condition P(z max , t) = for all t and the initial 
condition P(z, 0) = S(z) for z < z max . 

The presence of the sink at z = z max explains a phe- 
nomenon barely noticeable in Fig. 7, but exemplified 
in Fig. O namely the departure of the observed average 
(In T) from the usual DMPK result (InT) = -t. This de- 
parture is due to the fact that the numerically computed 
average only takes into account the stationary solutions 
which -as will be seen from Eq. (fTT))) - become less and 
less numerous when t increases. Hence what is computed 
numerically is the average of InT = —z over the distri- 
bution (TT6"1) . This reads 



zP(z,t) dz 



1 + erf 



- t 



Isft 



erfc 



t ~\~ z ma 
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where the error function is defined by 



erf (a;) 



2 f x 

7= / cxp(-y 2 )dy 
/7r Jo 



(77) 



(78) 



and erfc(a;) = 1 — erf(x). 

Expression ([77|) is compared in Fig. [9] with the results 
of a numerical simulation performed in the case V/c — 
450 (corresponding to z max = 9.43) for 10000 random 
potentials Us of type (|21[) characterized by ng£ — 0.5 and 
£/b= \/2 (leading to Li oc (k) = 100^). The agreement is 
seen to be very good. Since the sink cuts the solutions 
which are strongly scattered by the random potential, the 
remaining stationary states have a higher transmission 
coefficient. This effect increases with the sample length 
L, which explains the behavior of the curve in Fig. [9] 

As an other test of the validity of our approach (which 
amounts to model the upper boundary z max by a perfect 
sink and to neglect nonlinear deformations of the density 
pattern of the flow close to the threshold) we now deter- 
mine the probability for a trajectory to remain below the 
boundary. This survival probability reads 



P{z,t) dz 



2z n 



t) with a suitable prefactor. This yields the dis- 
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■ erfc 
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FIG. 9: (Color online) (InT) = — (z) plotted as a function of 
t — L/Li oc (k) in the case of random potentials Us character- 
ized by n s £ = 0.5 and £/b = y/2. The curve is drawn in the 
case V/c = 450. The red solid line is the numerical result and 
the black dashed line is the analytical result (|77|l . The straight 
(thin dashed) line is the usual DMPK result (In T) = -t [Eq. 

CD)]- 



As anticipated, P s (t) clearly decreases from 1 (at t = 
0) to (for large t). The knowledge of P s (t) allows to 
determine the value L* of the length of the disordered 
region beyond which most of the random realizations lead 
to a non-stationary flow of the condensate. We can, most 
conveniently, define L* through the condition 



P s (t*) = 1/2 



(80) 



with t* = L*/L\ oc . This leads to the implicit equation 
for the threshold value t* : 



rf 



- t* 



*erfc 



+ t* 



(81) 



This equation can be explicitly solved in the limiting case 
of large z max . As it is natural to assume that t* ought to 
be of the order of z maX ) which is the only relevant scale 
in this equation, we make the ansatz 

t* = Zmax + St (82) 

and assume (which is to be verified a posteriori) that 6t 
is of the order of unity, whereas z max 3> 1. This yields 
to lowest non vanishing order 
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t* 



St 
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FIG. 10: (Color online) Fraction of stationary trajectories 
P s (t) plotted as a function of the length L of the disordered 
region. The condensate flow through disorder potentials Us 
of type (|21[) was numerically computed for this purpose (red 
solid line), at parameters for which V 2 /c 2 = 2 x 10 5 . The 
black dashed line shows the analytical prediction of this sur- 
vival probability P s (L/L\ oc ) according to Eq. (|79[) . The verti- 
cal dashed line marks prediction (|86[1 for the threshold length 
L* at which P 3 (L* / ' L\ oc ) = 1/2 (horizontal dashed line), 
namely L* / L\ oc = 8.433 . 



for the right-hand side of Eq. (fSTj) . This finally results in 

6t = -l + 0(z~l x ) . (85) 

Neglecting terms of the order of z^^, we therefore obtain 
for the threshold length 



DL 



loc 



^loc 







Kl6c 2 ) 


- 1 





(86) 



We emphasize that Eq. (|86|) holds for z max ^> 1, i.e., 
for ln(F 2 /16c 2 ) > 1 [see Eq. 1(75)1] ■ This is much more 
restrictive than the condition V c which is assumed to 
hold true when deriving Eqs. (f7rJ|) and ((HJ). 

Figure [TU1 shows a comparison of the analytical predic- 
tions ((79)) and ((86)) with numerical data obtained from 
the integration of the time-dependent Gross-Pitaevskii 
equation ((1]). The condensate flows through a disorder 
potential U s of type j2T)) with V 2 /c 2 = 2 x 10 5 . We see 
that the fraction of stationary trajectories P s (t) is very 
well described by Eq. (fTTJ)) . and that the approximate ex- 
pression ((56")) predicts very well the length L* at which 
the crossover length from stationary to time-dependent 
flow occurs. 

For velocities V not extremely large compared to the 
speed of sound, the condition z max ^> 1 will not be ful- 
filled and estimate ((5fJ)) will not be valid, while A max ^> 1 
might still hold and the average evolution of the sys- 
tem might still be fairly well described by the simplified 
Fokker-Planck equation ((75)) . In that case, the implicit 
equation (|8ip has to be solved numerically. In Figure 
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FIG. 11: (Color online) L* as a function of V in dimensionless 
units. The black dashed line corresponds to the solution of 
Eq. (|8ip and the red solid line corresponds to the value of L* 
extracted from numerical simulations for a potential Us with 
the same characteristics as in Fig. [T](see the text). 



Q] one can see that the numerical solution of Eq. (jSTj) 
(yellow solid line) provides a very reasonable estimate of 
the boundary between the bright supersonic region (sta- 
tionary flows) and the dark time-dependent region, in a 
regime of not extremely large V/c, where Eq. ([55)1 fails 
to properly predict the threshold length L* . The simula- 
tions are performed by solving Eq. |T]) numerically using 
a potential of type Us (characterized by (Ug)/fi — 0.025 
and n 5 £ = 0.5). For each V and L we consider 100 re- 
alizations of such a potential and statistically determine 
the quantity P s , i.e., the fraction of stationary solutions. 
P s is plotted in Fig. Q] using a gray scale (dark, P s = 
; light blue/gray, P s — 1) as a function of the normal- 
ized variables L/£ and V/c (this normalization rescales 
interaction effects) . The qualitative agreement of Figure 
[T] is made quantitative in Figure [TTJ In this figure the 
numerical solution of Eq. (|8ip is compared with its de- 
termination extracted from numerical simulations in the 
supersonic regime. More precisely, the solid (red) line in 
Fig. [TTlis simply the contour P s = 1/2 in Fig. [TJ This 
corresponds exactly to definition (|80[> of L*. The agree- 
ment between the numerical result and the theory of the 
present section [dashed curve, solution of Eq. ([51])] is 
seen to be excellent [Hf| . 

We conclude this section by emphasizing that the exis- 
tence of an upper threshold L* corresponding to lengths 
of the disordered region beyond which most of the flows 
are time-dependent is a genuine nonlinear effect [absent if 
one sets g = in Eq. (JTJ)] . Actually, whereas interactions 
only weakly modify the precise value of the localization 
length, the existence of the threshold L* is a remarkable 
qualitative effect induced by nonlincarity. Moreover, as 
illustrated in Figures [5] and [TUl this effect persists even 
in the limit V 3> c where naively one would expect no 
noticeable consequence of interaction. 



VI. EXPERIMENTAL CONSIDERATIONS 

On the basis of the results obtained in the previous 
section we present here what are the more favorable ex- 
perimental configurations for observing Anderson local- 
ization in an interacting Bose-Einstein beam. We also 
discuss a possible experimental signature of localization. 



A. Appropriate configurations for observing 
Anderson localization 

In the non-interacting regime the only condition for ob- 
serving Anderson localization in ID is that the size of the 
disordered region should be larger than the localization 
length. Then one can observe an exponential decay of 
transmission with a log-normal distribution (in the limit 
L > Lioc). 

The situation is more complex when interactions are 
turned on. What is particularly interesting is the inter- 
play between localization and superfluidity. Indeed these 
two phenomena are conflicting one with the other: su- 
perfluidity is the (counterintuitive) ability to pass over an 
obstacle without reflection whereas Anderson localization 
corresponds to a large reflection in a domain where one 
would expect almost perfect transmission. As a result 
of the interplay between these two extreme phenomena, 
and depending on the fluid velocity and on the sample 
size, the flow may be stationary and superfluid, dissipa- 
tive and time dependent or stationary supersonic (and 
also dissipative) [14j . Anderson localization does not oc- 
cur in the superfluid region (where the transmission is 
perfect) and either does not exist or cannot be clearly 
identified in the time dependent regime (where interfer- 
ence effects are washed out [13j): but is truly observed 
in the supersonic stationary regime, as demonstrated in 
Sec. B 

In that regime, a first experimentally relevant effect is 
the modification of the localization length with respect 
to its value in the absence of interactions. This effect is 
very well described by renormalizing the wave vector k to 
k [Eq. (|48|) ]. which means that part of the kinetic energy 
available to the flow is taken by interactions. However, 
as already discussed in Ref. [T3| this effect is only sizable 
in a regime where V is not too large compared to c, and 
is thus relevant only in the perturbative regime (cf. Fig. 

A second experimentally observable effect is the modi- 
fication of the localization length due to the correlations 
of the disordered potential. This is described by formula 
([48]) where C is the Fourier transform of the two-point 
correlation function of the disorder. For the different po- 
tentials considered here, (7=1 for a potential Us or is 
alternatively given by Eqs. (|34p. (pJB")) and (|4"Tj) for corre- 
lated potentials. Explicitly this yields 

L 1oc (k) = — , (87) 
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for a potential Us of type (|2Tj) ; 

k 2 

iioc(«) = — exp{2 K 2 ^ 2 } , (88) 

for a potential J7 G of type l[25]): 

k 2 

£ioc(k) = — exp{2/cl c } , (89) 

(7 

for a potential C/x, of type (f2"5"|) ; and 



for a potential J7 S of type (|3G[) (when < 1). The 
validity of these expressions has been tested in Sec. IV Bl 
In the non-interacting case (i.e., k = k), expressions ([57)1 
to ([201) correspond to a high energy limit and can be ob- 
tained through a first order Born expansion within the 
phase formalism of Refs. 0, [47]]- In all three cases, one 
sees that the localization length is drastically enhanced 
due to the non-zero correlation length with respect to the 
uncorrelated disorder, Eq. ([57)) . In the Gaussian and the 
Lorentzian cases the localization length scales exponen- 
tially with (k£ c ) 2 and k£ c , respectively [see Eqs. (|88I89|) ]. 
In the case of a speckle potential, the effect is even 
stronger: one sometimes speaks of an "effective mobil- 
ity edge" [I?], Hl| , meaning that beyond a critical wave- 
vector (or a critical velocity) the localization length ([9TJ)) 
is infinite. This is an artifact of the Antsygina-Pastur- 
Slyusarev formula P5|) which can be corrected by going 
to higher orders (see Refs. [1^, [MI ) : the corrections 
to this result give a localization length which is finite, 
but typically larger than any other relevant scale in ex- 
perimental systems. 

Hence, in all the cases the dependence of the localiza- 
tion length with respect k (i.e., with velocity) is amplified 
by correlations. Mathematically this is due to the fact 
that the denominator in the Antsygina-Pastur-Slyusarev 
formula (|48p for the localization length in presence of 
correlations tends to zero when k£ c 1. In order to 
minimize this effect one needs to impose the following 
condition: 

k£ c <1 or V <y c = -^ = c|- . (91) 

m £ c l c 

In the r.h.s. of Eq. (|9~lj) we replaced k by mV/h because 
in practice condition (|9ip is verified in regimes where 
V > 3 c, i.e., when the approximation n ~ fc is sound. 
Note that this condition is arbitrary and is only superfi- 
cially analogous to the 3D Ioffc-Regel criterion [62j . The 
latter defines a true mobility edge that separates a metal- 
lic from a localized phase whereas Eq. (|9~lj) only requires 
that the localization length does not get too large. Un- 
derstood in this sense, the criterion (|9"Tj) is exactly equiva- 
lent to the definition of an "effective mobility edge" some- 
times used in the literature. 

In the absence of interactions it is always possible 
(at least theoretically) to define a system with a length 




V* v c 

FIG. 12: (Color online) Schematic phase diagram in arbitrary 
units. The (blue) dashed line is the localization length L\ oc 
and the solid line is the threshold length L* . V c = h/m£ c is 
the typical velocity beyond which it is almost impossible to 
observe Anderson localization in a realistic system (see text). 
The (blue) colored zone corresponds to the region where An- 
derson localization can be experimentally observed in pres- 
ence of interaction. 



L > L\ oc which verifies (|91|) ; i-e., a system where one 
can observe Anderson localization. If we now turn on 
interactions, a major effect is the appearance of a length 
scale L* which signals the onset, for L > L* , of a regime 
of time-dependent flows (cf. Sec. IV Cp . In this regime, 
Anderson localization disappears, and the time-averaged 
transmission coefficient scales as 1/L [l3j]. This is the 
most spectacular effect of interactions in the transport 
properties of the system. In order to observe Ander- 
son localization, the system size should therefore satisfy 
L\ oc < L < L* . In practice, one should be in a regime 
of parameters such as illustrated in Fig. [12] the crossing 
L* > L oc has to occur at a velocity lower than V c . 

Based on the numerical solution of Eq. ([8~Tj) one can 
show that the crossing L* > L\ oc occurs at a velocity 
V* ~ 7.95 c [H[ (see also Fig. [T]). This condition only 
allows the system to reach a (stationary) regime where 
L = L\ oc . But if one wants to observe Anderson lo- 
calization one should be able to reach a regime where 
L* ^ L > 2 L\ oc say, in order to get as close as possible 
to the domain of log-normal distribution of transmissions 
still remaining in the region of stationary flows. This im- 
poses V/c > 20. This must be supplemented by condi- 
tion (|9ip. i.e. V/c < £,/£ c - Hence the correlation length 
£ c should be smaller or equal to £/20. Fig. [TBI shows the 
phase diagrams of a one dimensional interacting beam 
of condensed atoms moving through a speckle potential 
in this regime. For plotting this diagram one has gener- 
ated 16 random potentials and studied in each case if a 
stationary solution exists or not. The dark blue region 
corresponds to a domain where no stationary solution 
exists while the light blue one corresponds to a domain 
where all the potentials admit a stationary solution (the 
color code is the same as in Fig. [T]and is explained in Sec. 
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FIG. 13: (Color online) Phase diagram displaying the fraction 
of stationary trajectory P 3 for a beam with velocity V moving 
in a speckle disorder of extent L. The figure has been drawn 
for a potential 17s of type (|36|) characterized by £ c /£ = 0.05 
and a — 3.93 /x 2 £. The light blue region corresponds to a 
domain of stationary flow (P s = 1: 100 % of the solutions are 
stationary, see the explanation in the text); the dark region 
corresponds to time dependent flow. The curves indicating 
the values of L\ oc and L* correspond to the analytical results 
(1901) and to the numerical solution of Eq. (|81|) . 



IV C| . The region between L\ oc and L* in Figure [TBI is the 
region where one can observe Anderson localization. 

Let us now evaluate the orders of magnitude of the 
different parameters allowing one to reach the appropri- 
ate regime. For concreteness we consider a beam of 87 Rb 
atoms such as the one of the Atom Optics Group at Lab- 
oratoire Charles Fabry de l'lnstitut d'Optique. For a 
correlation length of 0.26 (im the velocity cut-off for ob- 
serving Anderson localization is roughly V c ~ 2.7 mm/s. 
Note that in Ref. [2| the velocity of the expanding con- 
densate is about 1.6 mm/s, i.e., smaller that V c as it 
should. If we use the parameters of Ref. [H| (ani ~ 0.25, 

V ~ 9 mm/s, a — 5.3 nm, n\ = 45 atoms/^m, c ~ 0.9 
mm/s and £ ~ 0.8 fim) it is impossible to satisfy con- 
dition (f9"Tj) because V ~ 9 mm/s > V c and also because 
£ c /£ ~ 0.3. However the Atom Optics Group has re- 
cently improved the sensitivity of its detectors which is 
now close to being able to detect a density as small as 1 
atom//im. This allows one to work with a smaller density 
and to improve the ratio t c /^ which can be tuned down 
to the value of 0.05. Then the localization length can be 
selected by tuning the speckle amplitude. For instance 
Lioc = 0.25 mm can be obtained by choosing (U s ) = 34 
Hz at V — 1.6 mm/s. These parameters are close to 
those used in Fig. [13] and are reachable experimentally. 
However, for observing Anderson localization one needs 
to keep the beam stable for almost 1 s (0.31 s if we want 
L = 2L\ oc ) whereas in the current experiment this is only 
possible during 0.1 s, hence it is still a matter of debate 
to decide if the observation of Anderson localization of 
a Bose condensed beam in the presence of interaction is 



within the reach of present time technology. 

It is also interesting to make a connection between the 
physics described here and the recent experiment obser- 
vation of Anderson localization of a condensate expand- 
ing in a disordered potential performed in the same group 
0. Contrarily to the propagation of a beam studied in 
the present work, Ref. 0] considers the spreading of a 
wave packet (initially at rest) in a speckle potential. Af- 
ter a first stage of expansion, mainly driven by interac- 
tions, the experimental cloud expands with a constant 
velocity V ~ 1.6 mm/s but the particle density and the 
sound velocity are functions of the position. Therefore 
it is not possible to place this experiment on a single 
point of the phase diagram displayed in Fig. [T2J How- 
ever one can evaluate the ratio V/c at different positions, 
that is at different L. For instance, at the typical value 
L = L\ oc , Fig. 2 of reference [2J allows to calculate the 
sound velocity as well as the healing length £, yielding 
the typical experimental values V/c ~ 12 and L/£ ~ 55. 
Moreover the ratio between the typical disorder ampli- 
tude (U s ) and the chemical potential /x = g no is the 
same in Fig. 1131 and in experiment: (U s ) / fj, = 5 (note 
however that in the experimental case this is the value 
of the local chemical potential that matters). Hence, al- 
though the experimental setup of Ref. 0] forbids a direct 
comparison with the results of the present work, the es- 
timates of the typical values V/c ~ 12 and L/£ ~ 55 
indeed locate the experimental system within the regime 
of Anderson localization of Fig. IT51 



B. Experimental signature 

Once the appropriate regime of parameters for observ- 
ing Anderson localization in a Bose condensed beam has 
been determined, it is also important to identify possible 
experimental signatures. In our theoretical approach we 
use the transmission coefficient T of the beam over the 
disordered region as the relevant parameter. However the 
measure of T might be experimentally involved, and we 
propose here to use an other related quantity, namely the 
rate of energy dissipation (2?]]: E = —VF& where 



F d = ( dxn(x,t) 

Jwl ox 



(92) 



is the drag force exerted by the beam on the obstacle. 
Definition ((92)) is quite natural: the force exerted on the 
obstacle is the mean value of the operator d x U over the 
condensate wave function. It is rigorously justified by 
the analysis of Ref. 



65] 



in terms of stress tensor. In the 
stationary case, changing integration from x to X in 
a simple integration by parts yields 



A/a 



(93) 



In the r.h.s. of Eq. (|9"3"]) we made use of relations |5]| and 
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It has been shown in section IIII Bl that E^ = and 
that under assumption (|11[) (small nonlinearity and arbi- 
trary transmission or weak transmission and arbitrary 
nonlinearity) one has (see Eqs. (|20[) and (fH))) ) E^ — 
2h 2 n 2 /m(R/T) which yields 

2H 2 k 2 R , x 

F d = n - . 94 

m 1 

In the regime R <C 1, using Eq. (|44[) one recovers from 
the perturbative result already obtained in Ref. [65| : 
F d = -2n m\U(2 K )\ 2 [ /H 2 0. 

The physics embodied in Eq. (|94|) is rather simple and 
it is worth spending some time to discuss it. Consider 
an incident beam of particles with density ni nc and mo- 
mentum p — —tin moving from +oo towards an obstacle 
at rest. Part of the particles is transmitted (a fraction 
T) and the other part is reflected (a fraction R). The 
collisions are elastic and each of the reflected particles 
experiences an exchange of momentum bp = 2Hk with 
the obstacle. During a time St there are iV co n collisions 
and by the law of action and reaction the obstacle expe- 
riences a force 

dp 2h 2 n 2 , 

F d = -JV col] -£ = n inc R . (95) 

or m 

In the r.h.s. of (|95|) . one has written that N co \\/5t is the 
flux of particles colliding with the obstacle i.e., ^ Rn lnc . 
Eqs. (|94|) and (|95|) are identical because what we call no 
is the downstream density of the beam (cf. Fig. [2]), i.e., 
precisely Tni nc - Depending on which quantity is held 
constant (no as in the present paper, or n; nc ) Eq. (|94[) or 
(|95| is more appropriate (cf. the discussion of the fixed 
input and fixed output problem in Ref. t 32]). This is 
somewhat reminiscent of the controversy on the correct- 
ness of the Landauer formula (see, e.g., the discussion in 
Ref. [13). 

On the basis of (|94|) one sees that the measure of E 
gives direct informations on the transmission of the in- 
teracting beam through the disordered region, allowing 
one to reveal in which configuration is the system. For 
instance in the localized regime the energy dissipation is 
high (cx 1/T) and grows exponentially with the size L of 
the disordered region, whereas in the perturbative regime 
E is much lower and scales as L. 



assumption (jlip . which holds in the following regimes : 
(i) small nonlinearity and arbitrary transmission or (ii) 
weak transmission and arbitrary nonlinearity. In other 
cases there is no obvious way to define the transmission 
of the non-linear beam because one cannot separate in 
the up-stream region an incident flow from a reflected 
one. However, our analysis in terms of E c \ and A (defined 
in Sec. IIII Bp is always valid, even when condition (fTTj) 
is not fulfilled. This just means that, out of regime (jlip . 
the connection (f20|) between A and T is invalid. But for 
instance this does not invalidate at all the analysis lead- 
ing to the DMPK equation ([55]) . and the experimental 
signature proposed in Sec. IVI Bl also remains valid even 
when it is not possible to properly define T. 

We note that the validity of the DMPK approach for 
non-interacting particles is a well established fact in the 
theory of disordered systems. What is achieved in the 
present work is its extension to the case of interacting 
particles. Other studies of Anderson localization in the 
presence of interactions have concentrated on the long 
time behavior of the time evolution of initial wave pack- 
ets [68] . Although those results are still a matter of active 
debate in the community, the results of the present work 
produce strong evidence of the existence of Anderson lo- 
calization for weakly interacting Bose particles (with ef- 
fective repulsive interaction) propagating through disor- 
dered samples of finite size L < L* . 

Although the present study leads to the important con- 
clusion that Anderson localization in the presence of in- 
teraction is possible, it is rather disappointing to remark 
that it can be clearly identified only when V > 20c, i.e., 
in a regime where interactions do not play a major role 
(see the discussion of Sec. IVI Al and also Ref. [3|). In 
this respect, the more interesting and new effect of inter- 
actions is the existence of an upper threshold L* for the 
length of the disordered region : for L > L* no station- 
ary flow is possible. As shown in Sec. IV CI L* is directly 
connected to the probability distribution of the param- 
eter A and to the localization properties of the system. 
It would be very interesting to lead a systematic study 
of the transmission in the interaction-induced time de- 
pendent regime (L > L*) where the numerical results of 
Ref. [HI indicate a power law decay of the time-averaged 
transmission, a signature generally considered of loss of 
phase coherence and onset of Ohmic behavior [6^, FtoI ] . 
Work in this direction is in progress. 



VII. CONCLUSION 

In the present work we have presented an analysis of 
the transmission of a weakly interacting Bose gas incident 
on a disordered potential. We have shown on the basis 
of numerical and analytical results that there is a regime 
of Anderson localization in this system and proposed ex- 
perimental signature of this phenomenon. In order to 
properly identify a "localized regime" we have studied 
the transmission coefficient and its probability distribu- 
tion. The transmission coefficient T is well defined under 
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APPENDIX A: DISTRIBUTION OF 
REFLECTION COEFFICIENTS IN THE 
PERTURB ATIVE CASE FOR A POTENTIAL OF 
TYPE (p5j) 

We give here a demonstration of the perturbative re- 
sults (|46|) , (|47|) and {48]) in the special case of a Gaussian 
disordered potential U g of type (f25|) . A simple way to 
obtain this results starts by noticing that any Gaussian 
noise verifying (rj(x)) = and (rj(x)r/(x')) — S(x — x') 
[and here we are specifically interested in r](x) that ap- 
pears in Eq. (|2"5|) ] can be written as (see, e.g., Ref. ]T\\ ) 



T)(X) 



lim ■ 

v — >oo 



(Al) 



j=-oo 



where Xj's are random positions uniformly distributed 
on the real axis with density v and mean spacing 1/v 
and tj = ±1 is a random variable (with (ej) = and 
(ei€j) = Sij). 

In order to calculate the probability distribution of the 
reflection coefficient R [whose value is given by (|44[) ] one 
should first consider the distribution of 



U g {2n) = dxU g (x)e 



h 2 ^ w(2k) Z 



lim _ 

u—>oo m Jv — ' 
v j=o 



(A2) 



The quantity U g (2n) as given by (|A2|) is formally equiv- 
alent to the position z of a particle performing a random 
walk in the complex plane after N — vh iterations. The 
particle is initially at the origin and performs jumps of 
constant amplitude s = ^-|ui(2k)| \fo~fv with random di- 
rection. Denoting by d 2 P = p(z, N)dxdy the probability 
to find the particle in the domain dxdy around z after 
./V steps, if N ^> 1 [which is ensured by taking the limit 
v — ► oo in (|A2[) ] , the central limit theorem yields 



p(z,N) 



1 



ttNs 2 



exp 



Ns 2 



(A3) 



It is then a simple exercise to get the distribution of |z| 2 . 
One obtains 



P{\z\ 2 1 N) = ^rexp 

(pr) 



\z\< 



<M 2 > 



(A4) 



where 



z\ z ) = s N = a(h /m) \w(2k)\ L . (A5) 



From relation (|44[) , Eq. (|A5p immediately yields the an- 
nounced probability distribution l|46p with 



2 > = ^R2 K )| 2 



(A6) 



For a potential U„ of type 



.7,1 2 



C [see Eq. {30])] 
and Eq. (IA6|) demonstrates in this case the validity of 
Eqs. {T7]) and 



APPENDIX B: DERIVATION OF THE DMPK 
EQUATION (f66|) 

In this appendix, we explain how to obtain the DMPK 
equation (|66[) starting from the discrete Langevin equa- 
tion ([621) . 

Let us consider a generic situation where A„ obeys a 
stochastic recursion relation of the type 

Ki+l — Ki = F(X n ,(n) ; (Bl) 

with uncorrelated random variables Cn : 

It is clear that under assumption (163]) , Eq. ([62]) is of type 
(|B1[) with all the odd N averages in (|B2[) being zero and 
C* 2 = 1/2 [cf. Eq. {Bl])]. 

Let P(A, n)d\ be the probability that X n lies in the 
interval A, A + dA. One can express P(X,n) as 



P(X,n) = (S(X n - X)) = 

This yields 

P(A,n + l)-P(A,n) 



dk ; 



ifc(A„-A) 



2tt 



(B3) 



2?r V 



dk 



E^M^i^SiXn-X)) 



(B4) 



Using the fact that A„ depends on the variables Ci, C2 
. . . Cn-i but not on ( n (as can be seen directly from (|Blj) ) 
one can write the last of Eqs. (|B4j) as 

P{X,n + 1) - P(A,n) = 



E^rrayl^'Cn))^,-)} ' (B5) 



In the case of Eq. {62]) one has F(X, C) = (1 + 2A) /k 2 6 2 + 
2 (A 2 + A) 1 / 2 ^/ nb and the successive moments of F read 



<^(A,Cn)> 



1 + 2A 

K 2 b 2 



(B6) 
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(F 2 (\,Cn)) = 



2(A 2 + A) 
K 2 b 2 



O 



i% 4 



(B7) 



with all the other moments being of order l/(re 3 6 3 ) or 
more, i.e., negligible in regime ([60)) . Eq. (|B5[) thus reads 



In the continuous limit, defining t = n/(n 2 b 2 ), the l.h.s. 
of Eq. (|B8[) is simply the first derivative of P with respect 
to t and (|B8|) reduces to Eq. (|66[) of the main text. 



K 2 b 2 [P(A,n + l)-P(A,n)] 



2A)P] + ^[(A 2 + A)P] = 



A(A + 1) 



dP~ 



(B8) 
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